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VOLUME AND SURFACE INSTABILITY 
IN SLIDING PLASMAS 

by 

S. P. Talwar 

Goddard Space Flight Center 

SUMMARY 

The electrostatic instability arising due to a tangential 
discontinuity of velocity in two semi- infinite, homogeneous 
regions of a cold and collisionless plasma subject to a uni- 
form magnetic field is investigated by use of two fluid equa- 
tions. A general dispersion relation is derived allowing for 
interstreaming of charged particles in each region in addi- 
tion to a slippage of one region over the other. Special 
cases are then discussed and the criteria for stability, 
monotonic instability and growing wave instability are 
derived. 
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VOLUME AND SURFACE INSTABILITY 
IN SLIDING PLASMAS 


by 

S. P. Talwar* 

Goddard Space Flight Center 


INTRODUCTION 

Plasmas characterized by relative tangential motion between various layers are a common 
occurrence in many astrophysical situations, e.g., solar-wind magnetospheric boundary and 
coronal streamers moving through the solar wind. This had led some workers (References 1, 2, 
3, and 4) to investigate, in recent years, the instability which arises due to a tangential discontin- 
uity in velocities in hydromagnetic plasmas. These investigations made use of the collision- 
dominated single fluid hydromagnetic equations, although in a recent investigation (Reference 5) 
the corresponding problem for a plasma of weak collisions characterized by an anisotropic 
plasma pressure is worked out using the Chew, Goldberger, and Low equations (Reference 6). 

Quite often (e.g., in enhanced emission of charged particles at the time of a solar flare) situa- 
tions arise where a plasma, carrying a net volume current, penetrates another plasma thus leading 
to a configuration where one expects to find the joint effect of two- stream instability and Kelvin- 
Helmholtz instability. In the present paper we intend to investigate the instability of such a plasma 
configuration while neglecting the effects of finite temperature and collisions. 


EQUATIONS OF THE PROBLEM 

Consider two semi-infinite, homogeneous, electrically neutral plasmas slipping past one another 
at the common interface z = 0. Each plasma is regarded as cold and collisionless and subject to 
a uniform magnetic field parallel to the interface. The particles in each plasma region (z < o and 
z > o) are assumed to move parallel to each other in the initial state and at constant speeds, V 1 . and 
U 2 , respectively along the x -direction. Here j stands for either electrons or ions. The electrons 
and ions may move together in each plasma in which case there will be no initial current flowing 
in either plasma, and the configuration of sliding plasmas will correspond to the classical Kelvin- 
Helmholtz problem as applied to ionized gases. The electrons and ions may, however, be taken to 
move at different speeds (but each spatially constant) in each plasma region. In this case there 
will be present a contra- streaming of charged particles in the body of each plasma in addition to a 
tangential discontinuity of velocities for each charged species at the interface between the two 
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plasmas. Here, the mean electric current is not zero in the initial state, and one should expect a 
simultaneous presence of both the contrastreaming instability (the volume instability) and the Kelvin 
instability (the surface instability). We shall neglect the self magnetic field arising due to the in- 
itial current as compared to the prevailing magnetic field of external origin. The prevailing mag- 
netic field may, in general, be inclined at a certain non- zero angle (but still parallel to the inter- 
face) to the direction (x-axis) of free streaming of charged particles in the initial state. If the 
magnetic field is assumed along the direction of free streaming there will be no initial electric 
field in either medium. On the other hand, if the prevailing magnetic field is, in particular, taken 
normal to the free streaming, then each plasma region will be characterized by an initial constant 
electric field in a direction normal to both the streaming motion and the magnetic field, i.e., nor- 
mal to the interface (the z -direction). In case the ions and the electrons are not moving together 
in such a configuration, we may reconcile this by arguing that electrons (streaming through sta- 
tionary ions) acquire a steady drift velocity in the crossed e and B fields and that the ions do not 
partake in this motion. Such a situation may arise in various physical situations (e.g., in the ion- 
osphere in the form of the "electrojet" where a strong Hall current flows perpendicular to the 
earth's magnetic field, and in PIG type discharges (Reference 7) in the laboratory). 

To investigate the stability of the above mentioned configuration, we will develop the perturba- 
tion equations, for the initial state which is defined as follows: 

B 0 = 


= 

Po = 

For a zero temperature plasma, in which collisions are disregarded, the two fluid perturbation 
equations are written as, 


[Box' B 0y , 0] 

Kj ■ o, o] 


constant 


(i) 


(A + u oi • v ) 


e Uj X B 0 U X b 

in” V* + — + ~ 


( 2 ) 
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V X E 


(6) 


and 


_ IjL Z 

c <9t b 


V • E 


477 L e i n i 


(7) 


Here u. and n j denote, respectively the perturbation in velocity at a point and the number density 
of the j th class of particles (electron and ion). The quantities b, andT denote the perturbation in 
the magnetic field and current density at a point, E stands for the electric field at a point in the 
medium, and the equilibrium quantities are shown with a suffix ‘0’ . For simplicity we shall, in 
what follows, restrict our discussion to electrostatic perturbations only so that we may regard the 
perturbation b in magnetic field to be zero. Thus the last term in Equation 2 vanishes and Equa- 
tion 6 is written as, 


so that 


V x E = 0 , 


E = -V0 


(8) 


(9) 


By coupling Equations 7 and 9, we get 


V 2 0 + 4 rr 0 . (10) 

If we assume the components of perturbation to vary with x, v, z, and t as (some function of z) 

• exp i (k x x + k y y - at) , we can rewrite Equations 2 and 3 as 


^-ioj + ik * U Qj ) n. - “N Q V • ^ (11) 

and 

+ ik 1 U oj ) u. +0. xu. “ " m7 . (* 2 ) 

where O j stands for the Larmor frequency e. i/m. c for the charged particles (electrons, -e 

and ions, +e), and has components O jx and n. y which correspond, respectively, to the magnetic field 

components, B 0x andB 0y in x- and y -directions. 


3 



Using Equation 11 we rewrite the perturbation equation for the electric potential, 0, as 


V 2 \p - 4" N 



0 . 


( 13 ) 


We now have to evaluate an expression for V • IT in terms of 4> to be substituted in the above equa- 
tion so as to obtain the final equation in a single perturbation quantity «, p. 

Taking the divergence and the curl of Equation 12, and simplifying we obtain 


■ 9 ( r r; V 

-1 



/ _ _ „ (n. • v)(n. • Vi//) 

( ■ t -T u r?i , V J A J _ Y i 

n i r _k ' u oJ 

V ■ u. 
j 


m . 
j 

V / Y l-iu + ik • U oj J 


Eliminating V • u j from Equations 13 and 14 we obtain 



where co denotes the plasma frequency (477 n 0 e 2 /m.Y /2 for the j th class of particles. 

p j \ / 

With the spatial dependence of perturbations envisaged, we may rewrite Equation 15 as, 



where D stands for d/dz . 


The above equation with u oj = 0 and D 2 = -k z 2 gives the results for electrostatic plasma oscil- 
lations in a homogeneous static plasma subject to a uniform magnetic field. In particular it leads 
to the following well-known results. 


a 


,2 



(Longitudinal propagation) 


and 


1 + 


n . 2 


(Transverse propagation) . (17) 
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Again, we may deduce from Equation 16 the results of two- stream instability for a homogeneous 
plasma carrying a uniform magnetic field along the direction of streaming, (Reference 8) and 
transverse to the direction of streaming (Reference 9) by taking u oj ^ o and d = 0. 


BOUNDARY CONDITIONS AND DISPERSION RELATION 

For a configuration of two plasmas sliding past each other Equation 16 holds for both plasma 
regions (z < o) with respective values for the physical parameters involved in the equation. The 
respective solutions of Equation 16, bounded at z = ±°°, for the two plasma regions z < o and z > o 
are written as, 


and 


. _ m z 

'Pi (z) - Ae 


(z <0) 


> (18) 


P 2 (z) = Be m22 (z > 0) , 


J 


where m t and m 2 are chosen so as to have positive real parts, and are given by 

/k 2 


k 2 



„(C •»„)’> 


(" ■ E iJ .,J 2 [ n ,? - (“ - k • B, i) ’] 

- (c.-k - U ,,) 2 


(19) 


here i = 1, 2 for the two regions z <0 and z > 0 , respectively. 

We need to satisfy appropriate boundary conditions at the common interface between the two 
plasmas. The conditions are: 

(1) Equation 8 gives the boundary condition that ^ should be continuous across the interface. 
This yields 


A = B . 


( 20 ) 


(2) Equation 13, when integrated across the interface, yields the condition 


IVi 



(- 


477 N 
i + ik • 


01 



u . 
J z i 
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( 21 ) 
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The expression for u . z can be derived from Equation 12 written in component form. After some 
simplifications we obtain, 

(-io> + ik • U Q j) 2 + f2j 2 "j u j z = " [(-ico + ik • U oj ) + i(kxH.)^J • (22) 


Substituting Equation 22 into Equation 21 and making use of Equations 18 and 20 we obtain 



Here m 1 and m 2 are given by Equation 19 with necessary suffixes. Equation 23, with values for 
m 1 and m 2 substituted from Equation 19, constitutes the dispersion relation determining the stability 
of the configuration under investigation. It is rather unwieldy for discussion in the general case, 
and we shall therefore discuss some special cases taking the propagation vector k parallel to the 
streaming motion. For a propagation vector transverse to the direction of streaming, the config- 
uration is equivalent to a static state and is thus stable in the absence of any discontinuity in mag- 
netic field and plasma density at the interface. 


FIELD-FREE SLIDING PLASMAS 

If both plasma regions (z > o) are devoid of any magnetic fields, andm 2 are each equal to k 
and the dispersion relation, Equation 23, reduces to 

_I , _l^!_ , — ( 1 , -^'Ml , , 24) 

(“-«„)■ t»-ku,.)» » Vj 

where m and M, respectively, stand for electron and ion masses, andN x and n 2 are the particle 
densities in the two regions. In Equation 24 we have considered the propagation vector k along the 
direction of streaming (namely the x-axis). For simplicity we may further restrict ourselves to 
a single homogeneous plasma (n x = n 2 ) wherein the ions are at rest and the electrons stream past 
them with equal and opposite speeds in the two semi- infinite regions z < 0 and z > o. Thus, with 
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u ii = U 2i - 0 and U le - -U 2e 


u 0 , Equation 24 reduces to 


1 2 ^ [(^" ku o) 2 + ( a; + ku o) 2 + 2 Mcu2 ] ^ 

This is a cubic equation in co 2 which does not possess any negative real root for co 2 leading to mono- 
tonic instability in case the inequality k 2 U 0 2 >o; p 2 is satisfied. Equation 24 may be compared with 
the corresponding formula for the two- stream instability in a homogeneous current-carrying 
plasma in the absence of any magnetic field. Measurings, k u 0 in units of s we may rewrite the 
Equation 25 as 


F ° (X ' y) " 1 " [x 2 + (x-y) 2 ] 

where 


It may be noted that our Equation 26 differs in the last term on the left hand side from the conven- 
tional formula for two- stream instability in a homogeneous plasma devoid of any relative slippage 
between its two halves. The function, F 0 (x, y), and the corresponding function, F(x, y) without the 
last term 2xy/(x 2 - y 2 ) 2 , are plotted in Figures 1 and 2 for y = 0.5 and 1.5 respectively. The curve 
F 0 (x, y) , symmetric in x, is always less negative than the curve F(x, y). For y > l, the plot of both 
functions F and F 0 demonstrates that all roots (four for f and six for f 0 ) are real, leading to stable 
oscillatory modes. For y < 1, on the other hand, the function f corresponds to an over stable situa- 
tion (instability through growing waves) whereas the function f 0 leads to a monotonically unstable 
situation corresponding to a negative real quantity. Let us now assume that the ions and the 
electrons move together so that there is no initial current in either plasma region. In this case, 
the dispersion relation, Equation 24, reduces to, 


2xy 


(x 2 - y 2 ) 2 


= 0 


ku n 


and 


m_ 

M 


(26) 


(27) 


1 


H i + " l 

(c-kUj 2 (w~kU 2 ) 2 


(28) 


where 


+ Cl) ' 


There is no contrastreaming of particles in the body of either plasma, and the only discon- 
tinuity in motion is present at the interface. Equation 28 for co p ^ and \J 1 - -u 2 r U 0 
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Figure 1— Plot- of the functions F 0 (x, y) and F(x, y) versus x for y = 0.5. 


2 ' 



This formula is identical to the one obtained by Barston (Reference 10) and shows that a tangential 
discontinuity in a field-free homogeneous plasma is stable or unstable monotonically depending on 



whether the relative speed 2U 0 is more or less than the critical value defined by 

H>). = • (3°) 

It may be noted that the classical Kelvin- Helmholtz instability in hydrodynamics is recovered 
in the limit ^ Further, the Equation 28 shows that there is no possibility of growing wave 

instability in the configuration. In Figure 3 we have shown in curve a, the dependence of the root 
o) 2 /oo 2 (corresponding to negative sign in Equation 28) on the parameter ku 0 /co . For a given 
tangential speed the wave numbers of perturbation lying between the range 0 and k* (Equation 30) 
are monotonically unstable, there being a mode of maximum instability for an intermediate wave 
number defined by 


k 


m ax 



with maximum growth rate given by 



(31) 


(32) 


SLIDING PLASMAS WITH A UNIFORM FIELD ALONG STREAMING 


If we suppose that the prevailing magnetic field is parallel to both the direction of streaming 
of particles in both plasma regions z < o and the propagation vector k, then the terms involving 
(k x H ) disappear in Equation 22. The disper- 
sion relation, Equation 23, is thus written as, 



Substituting for m 1 and m 2 from Equation 18 we 
obtain the dispersion relation as, 




kU 0 /up e 

Figure 3 — Plot of the dependence of the root <l > 2 / co 2 
(from Equation 36) on the parameter kU 0 /aj . Curve a e , 
n y% = 0; curve b,n e /c.) pe = 0.5; curve c,n*/a) p =1.0; 
and curve d,f2 e /a> = 2.0. * 
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where 


"I,- - « - ku ij 


■\ 


and 


( 35 ) 


"2j - “ " kU 2j 


Electron Oscillations in a Single Plasma 

For convenience in handling the dispersion Equation 34 we may restrict our discussion to the 
case of electron oscillations alone and thus regard the ions to be unperturbed due to their heavy 
mass. In this case the summations in Equation 34 disappear and we get a simplified dispersion 
relation for a homogeneous plasma with a uniform magnetic field, so that the two halves of the 
plasma slide past each other along the field lines. This is written as, 


w 4 - ^2k 2 u 0 2 + n e 2 + to 2 ^ OJ 

Here we have again taken U x = -U 2 = U 0 
tion of Equation 36. 


2 + [ k4 U 0 4 - (fl e 2 + % 2 J(k 2 U 0 2 " n e 2 / 2 )] 

and cancelled a non- zero factor, (v 


2 
e 2 


o . 



(36) 

in the deriva- 


Equation 36 suggests that the configuration envisaged here does not show any overstability 
(growing wave instability corresponding to oj 2 a complex quantity) in case the inequality fi e 2 < 2k 2 u o 2 
is satisfied. The configuration is then stable or monotonically unstable depending on whether the 
inequality 


(vo 


2k 4 U 0 4 


(2k 2 U 0 2 " fi e 2 ) 


or 


> (37) 


2k 4 U 0 4 

( 2 k 2 u 0 2 -n e 2 ) 

is satisfied. The configuration will be overstable only if H e >w and the relative velocity satisfies 
the inequality 


8k 2 U 0 2 < (fl. 2 "« p a J • 


(38) 
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The dependence of one root, to 2 (in units of oj 2 ), likely to show instability on ku 0 (measured in 
units of to ) is exhibited in Figure 3 for n L = 0.5, 1.0 and 2.0 (curves b , c, and d respectively). 

P e e / P e 

Curve a represents the case of zero magnetic field. Clearly the configuration with field is unstable 
monotonically for a range of relative speeds, given by the following expressions, in case the elec- 
tron Larmor frequency is less than the plasma frequency. 


and 



The configuration is thus stable outside this range of relative speeds. The range of unstable rela- 
tive speeds shrinks to zero as the electron Larmor frequency approaches the plasma frequency. 
Thereafter the configuration is either stable or exhibits overstability for low enough relative 
speeds. For ft - 2 o> , for example, the configuration is overstable up to ku 0 - 0.6 a; and stable 
beyond. 

Thus we conclude by saying that the effect of magnetic field along the streaming direction is 
stabilizing in that the range of unstable tangential speeds is narrowed and the growth rate for any 
unstable wavelength is suppressed, till beyond a critical value of magnetic field (corresponding to 
: w ) the instability appears in the form of overstability for low tangential speeds only. For 
infinitely large magnetic fields, therefore, the configuration is unstable through over stability. The 
existence of instability for infinitely large magnetic fields was established earlier by Harrison 
and Stringer (Reference 11). Their analysis, however, did not apply to finite values of magnetic 
fields. 


Both Ions and Electrons Perturbed 


The assumption that only electrons are perturbed and not the ions is justified in case the dop- 
pler shifted frequencies, a> x and a > 2 , for either region are much larger than the ion cyclotron fre- 
quency. In general, ions also contribute and the general dispersion relation, Equation 33 is a 16th 
degree polynomial in o>, the parameter deciding the question of stability of the configuration. To 
reduce the degree of the dispersion relation when both ions and electrons are perturbed we may 
consider the following special situations. 
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No initial electric current in uniform plasma 


If the ions and electrons are taken to move together along the magnetic field lines, the disper- 
sion relation for the propagation vector along the streaming direction is written as, 


[— — 


— 


n 2 L 2 +n 2 ') 



L 2 +n. 2 ) 

e \ p i 1 / 


k<u 0 4 -fi e 2 ( 2 k 2 u 0 2 -n. 2 ) 

\ Pj i / 

k2U ° 2 + (a, 2 + fi 2 ) 

CD 2 + 

1 

L 2 + n 2 ) 

\ P e e / 


_ 

V Pe e l 


Here we have taken the doppler- shifted frequencies (a>-ku 0 ) and (o> + ku 0 ) for the two plasma regions 
to be very much less than the electron cyclotron frequency, . Thus, we find that in the approxima- 
tion mentioned above, the configuration of slow relative slippage in a uniform current free plasma 
is stable in the presence of a uniform magnetic field, along the streaming motion, and for all tan- 
gential speeds except for a narrow range of tangential speeds lying between the values correspond- 
ing to the positive and negative signs before the square root term in the following expression. 



It may be noted that there is no growing wave instability possible in this approximation. 


(41) 


Uniform plasma carrying an initial current 

Let us now suppose that in the two halves of a single plasma permeated with a homogeneous 
magnetic field, the electrons stream past the stationary ions with equal and opposite speed along 
the direction of the prevailing magnetic field. In this case the dispersion formula, Equation 34, 
for the case of slow relative velocities reduces to, 


2 co 2 (nf + k 2 U 0 2 ) 




+ 


2oo 2 (k 2 U 0 2 




0 ■ (42) 


This is a cubic in ^ 2 and so must possess at least one real root, which for monotonic instability is 
required to be negative. A necessary condition for monotonic instability to arise is that the elec- 
tron plasma frequency exceeds the electron Larmor frequency. In that case the configuration is 
unstable monotonically for a small range of relative speeds given by the expression 




( 43 ) 
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It may be noted that in the case of electrons streaming through stationary ions, the monotoni- 
cally unstable range of relative speeds (obtained from Equation 43) is very much less than the un- 
stable range of speeds when ions are regarded as unperturbed (Equation 39). 

The Equation 42 reduces to a quadratic equation in oP for frequencies much less than the ion 
Larmor frequency and in this limit we conclude that the configuration is monotonically unstable for 
relative speeds in the range defined by Equation 43 and is unstable through overstability if the fol- 
lowing inequality holds 


k 2 u 


2 

0 



(44) 


Clearly a necessary condition for growing wave instability to occur for low relative speeds is that 
the electron plasma frequency is less than the electron Larmor frequency. Thus, the configura- 
tion is monotonically unstable, though only for a narrow range of tangential speeds, if w > Q e ; and 
overstable for low enough tangential speeds if oj < . 

In Figure 4 we have shown the variation of one root (the less positive) of the quadratic equa- 
tion in w 2 , as obtained by neglecting oP in comparison to n . 2 in the term involving ^ p 2 in Equation 
42, with the parameter ku 0 for a few values of Q. e . The curves a, b, c are respectively for 
fi e /aj p - 0.5, 1, and 2. In the last value the configuration is overstable for relative speeds below 
a critical value defined by ku o - 0.6 & . 


SLIDING PLASMAS WITH TRANSVERSE 
MAGNETIC FIELD 

In case the two plasma regions (z > o) are 
assumed to carry a uniform magnetic field 
transverse to the direction of streaming in each 
region and the propagation vector k, and m l and 
m 2 (Equation 19) are each equal to k and the 
dispersion relation, Equation 23, is thus written 
as 

op 



("~ ku lj) [ n jl + ("~ ku lj)] 

.y r “ Pj l T = o (45) 

/ >-fcu 2i ) r°j2 - (“- ku 2i)i 



Figure 4 — Plot of the dependence of the root 2 /% 2 
(from Equation 42) on the parameter kU Q /£> p . Curve a, 
= 0.5; curve = 1.0; and curve c, 

= 2 . 0 . 
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For electron oscillations alone in a single plasma it further reduces to a quadratic in cP , given by 


[< + k 2 u 0 2 + (n. - ku 0 ) 2 ] + ku 0 (n. - ku 0 ) [ % 2 + ku 0 (n e -ku 0 )] = o . 


(46) 


Here we have taken Uj — U 2 - U 0 . 

It follows from the Equation 46 that the electron oscillations, in a uniform plasma with its two 
halves slipping past each other in a transverse uniform magnetic field, do not lead to growing wave 
instability. They do, however, show monotonic increase in amplitude in case the following condi- 
tions are simultaneously met with, i.e., 


kU A 




and 


> ( 47 ) 


ku 0 (ku 0 -n e ) , 


so that the configuration is stable for all electron densities in case ku 0 < ft e . 

Alternatively, we could argue that the electron oscillations lead to a monotonic instability 
only for wave numbers of perturbation greater than a certain critical value given by the expression, 


K)* 




(48) 


If one incorporates mutual streaming of particles also, leading to an initial electric current 
in each region, the dispersion relation, Equation 45, complicated for the general case, reduces to 
a cubic in oP for the case of electrons streaming through stationary ions in a single plasma with 
equal and opposite speeds in the two regions z <0. This is written as, 


1 + 


(<v-- 2 ) 


[2k 2 u 0 2 + (n - ku 0 ) 2 ] + k 2 u 0 2 (n - ku 0 ) : 


<0 2 lco 2 + k 2 U 2 - kU. n ) = 0 

p \ 0 0 e/ 


(49) 


This equation gives the following condition for monotonic instability in the medium. 


J 0 > ~2 



4 co 2 n . 2 

1/2] 


p e 1 


f 

1 + n 2 (n . 2 + oj 2 ^ 

e V 1 p ij 

1 

L 1 

—I 

J 


(50) 
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A comparison of the two equations, 48 and 50, shows that the critical wave number beyond which 
the configuration is unstable monotonically, is somewhat greater when perturbed motions of both 
electrons and ions are taken into account in a current carrying plasma than without it. 
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